Scaling at the chaos threshold in an interacting quantum dot 
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The chaotic mixing by random two-body interactions of many-electron Fock states in a confined 
geometry is investigated numerically and compared with analytical predictions. Two distinct regimes 
are found in the dependence of the inverse participation ratio in Fock space I on the dimensionless 
conductance of the quantum dot g and the excitation energy e. In both regimes 1^1, but only the 
small-p regime is described by the golden rule. The crossover region is characterized by a maximum 
in a scaling function that becomes more pronounced with increasing excitation energy. The scaling 
parameter that governs the transition is {e/g)\ng. 
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The highly excited atomic nucleus was the first exam- 
ple of a quantum chaotic system, although the interpre- 
tation of Wigner's distribution of level spacings as a 
signature of quantum chaos came many years later, from 
the study of electron billiards 1^. While the spectral 
statistics of the nucleus and the billiard is basically the 
same, the origin of the chaotic behavior is entirely differ- 
ent 1^ : In the billiard chaos appears in the single-particle 
spectrum as a result of boundary scattering, while in the 
nucleus chaos appears in the many-particle spectrum as 
a result of interactions. 

The study of the interaction-induced transition to 
chaos entered condensed matter physics with the real- 
ization that a semiconductor quantum dot could be seen 
as an artificial atom or compound nucleus A partic- 
ularly influential paper by Altshuler, Gefen, Kamenev, 
and Levitov ||^ studied the interaction-induced decay of 
a quasiparticle in a quantum dot and interpreted the 
broadening of the peaks in the single-particle density of 
states as a delocalization transition in Fock space. Differ- 
ent scenario's leading to a smooth rather than an abrupt 
transition from localized to extended states were consid- 
ered later [|[j8) . Recent computer simulations [^,0 also 
confirm the smooth crossover from localized to delocal- 
ized regime for quasiparticle decay. 

As emphasized by Altshuler et al. Q, the delocal- 
ized regime in the quasiparticle decay problem is not yet 
chaotic because the states do not extend uniformly over 
the Fock space. One may study the transition to chaos in 
the single-particle density of states, but theoretically it 
is easier to consider instead the mixing by interactions of 
arbitrary many-particle states. This was the approach 
taken in Rcfs. [p|J^jri|-[T^ , focusing on two quantities: 
The distribution of the energy level spacings and the in- 
verse participation ratio (IPR) of the wave functions in 



Fock space. Both quantities can serve as a signature for 
chaotic behavior, the spacing distribution by comparing 
with Wigner's distribution ^ and the IPR by compar- 
ing with the golden rule (according to which the IPR is 
the mean spacing 5 of the many-particle states divided 
by the mean decay rate F of a non-interacting many- 
particle state [l^]). Two fundamental questions in these 
investigations are: 1) What is the scaling parameter that 
governs the transition to chaos? 2) How sharp is the 
transition? 

In a recent paper p^ ] one of us presented analytical ar- 
guments for a singular threshold governed by the scaling 
parameter x = (e/gA) \ng, where A is the single-particle 
level spacing, e is the excitation energy, and g is the con- 
ductance in units of /h. (Both e/A and g are assumed 
to be 1.) In contrast, Georgeot and Shepelyansky |l^] 
argued for a smo oth crossover governed by the parame- 
ter y = (e/ff A)i,/£/A. (The same scaling parameter was 
used in Refs. pjl^].) The parameter y is the ratio of the 
strength V ^ IS./ g of the screened Coulomb interaction 
and the energy spacing A2 '-^ (e/A)~^/^A of states 
that are directly coupled by the two-body interaction Q . 
The parameter x follows if one considers contributions to 
the IPR that involve the effective interaction of 2, 3, 4, . . . 
particles. Subsequent terms in this series are smaller by 
a factor (ln5r/5)A„/A„+i, where A„ (e/A)-"+^/2A 
is the spacing of states that are coupled by an effective 
interaction of n particles [|4j. (The large logarithm \ng 
appears in the expansion parameter because of the large 
contribution from intermediate states whose energies are 
close to the states to be mixed.) 

The purpose of this paper is to investigate the 
interaction-induced transition to chaos by exact diago- 
nalization of a model Hamiltonian. We concentrate on 
the IPR because for that quantity an analytical predic- 
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tion exists for the e and f/-dependence. (There is no 
such prediction for the spacing distribution.) The numer- 
ical data is consistent with a chaos threshold at a value 
of X of order unity. 

The model for interacting spinless fermions that we 
study is the layer model introduced in Ref. and used 
for the quasiparticle decay problem in Ref. The 
Hamiltonian is H — Hq + Hi , with 

^o = 5Z^jS^S-' ^1 = H "^vMcjclc^c^. (1) 

j i<j,k<l 

The single-particle levels Sj are uniformly distributed in 
the interval [{j — 5) A, (j + ^)A]. The interaction ma- 
trix elements Vij^ki are zero unless i,j,k,l are four dis- 
tinct indices with i + j = k + I. The (real) non-zero 
matrix elements have a Gaussian distribution with zero 
mean and variance = (A/g)^. The Fock states are 
eigenstates of Hq, given by Slater determinants of the 
occupied levels fci, A:2, ^3, ■ • ■• The interaction mixes Fock 
states for which J2p equals a given integer. (Without 
this restriction the model is the same as the two-body 
random-interaction model introduced in nuclear physics 
1^,0.) The excitation energies of the states with given 
^1, ^2, fcs, • ■ ■ lie in a relatively narrow layer (width of or- 
der j^/'*A) around the mean excitation energy jA. The 
number of states in the j-th layer is the number of par- 
titions V{j) of j. For our largest j — 26 this number is 
7^(26) = 2436, which is still tractable for an exact diag- 
onalization. Without the decoupling of the entire Fock 
space into distinct layers, such large excitation energies 
would not be accessible numerically. The layer approx- 
imation becomes more reasonable for larger g, because 
then y ^ A so that states from different layers may be 
regarded as uncoupled. 

The inverse participation ratio / = X^mK*^!™)!^ 
the eigenstate \a) oi H is the inverse of the number of 
eigenstates \m) of TJq that have significant overlap with 
I a). We calculate / as a function of g for different layers 
J, corresponding to a mean excitation energy e = jA. 
The IPR fluctuates strongly from state to state and for 
different realisations of the random matrix H . We cal- 
culate the averages /, \/I and In/ where the overline 
~ indicates an average both over the V{j) states |a) in 
the j-th layer and over some 10'^ realisations of H . We 
first consider the logarithmic average In /, for which the 
fluctuations are smallest. 

In Fig. 1^ we have plotted the numerical data for the 
g-dependence of In/, for different values of e/A. In or- 
der to compare with the analytical prediction of Ref. 
[Q, we have rescaled the variables such that Fig. |l| be- 
comes a plot of —y~^h\I versus x. The prediction is 
that, in the thermodynamic limit the scaling func- 
tion F{x) = — y~^ln/ depends only on x for a; < 1. 
This scaling behavior can not be checked directly because 
finite-size effects introduce an additional e-dependence 
into the function F{x). This is why we can not directly 
test whether x or y is the correct scaling parameter. For- 



tunately, it is possible to include finite-size effects into 
the scaling function and test the theory in this way. 

Applying the method of Ref. for the calculation of 
In/ one finds that the function F{x) in the thermody- 
namic limit has the Taylor series 

oc 

F(x) = -y-ihrr-;^c„a;", (2) 

with corrections of order 1/lng. All coefficients c„ are 
positive. The scaling behavior (||) is expected to be uni- 
versal (valid for any model with random two-body inter- 
actions), but the coefficients c„ are model specific. The 
first two coefficients for the layer model are 

CO = ^^^^ = 1.53 , ci = |,/^co = 3.95. (3) 
V37r 25 V TT 

Finite-size effects introduce an e-dependence into the 
coefficients, through the e-dependence of /C„ = 
(A/A„)(A/e)"~-^/^. The series expansion of F{x) in 
terms of the /C„'s is 

F{x) ^ 4{V2 - 1)V^K2 + 36(2 - V2)K3X + 0{x^). (4) 

For e/A ^ 00 we have (2/7r)^273 = 

0.5198, K-i 6V6/257r = 0.1871 and we re- 

cover the thermodynamic limit (||). For the excita- 
tion energies e/A = 15,20,22,24,26 of the simula- 
tion, one has K2 = 0.419,0.436,0.439,0.444,0.447 and 
K3 = 0.0414,0.0536,0.0577,0.0615,0.0648. The result- 
ing small-x behavior of the scaling function is plotted in 
Fig. |l| (solid lines) and agrees quite well with the numer- 
ical data. 

Analytically, the scaling function F(x) is only known 
for a; <C 1. In the simulation, we observe a maximum 
of — y~^ln/ at x ~ 1. The maximum becomes more 
pronounced with increasing excitation energy. We argue 
that it is a signature of the transition to chaos, because 
beyond the maximum, for a; > 1, the IPR is observed to 
follow the golden-rule prediction (see discussion below) 

/goldcn-rulc = C[f/^V{j)]-^g^. (5) 

This golden-rule prediction is shown dashed in Fig. 
with the coefficient C « 0.51 as the single fit parameter. 
(The smallest e/A = 15 was left out of the fit.) Note that 
—y~^ In /goidcn-ruio has a maximum for an IPR of order 
unity, hence in the regime of localized states. In contrast, 
the maximum in —y~HnI occurs when the IPR is ^ 1, 
hence in the regime of extended states. We now discuss 
the small and large-x regimes in some more detail. 
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FIG. 1. Average logarithm of the inverse participation 
ratio / as a function of the dimensionless conductance g, in 
rescaled variables. The different sets of data points follow 
from the layer model for different excitation energies j = e/A. 
Statistical errors are smaller than the size of the markers. 
The straight solid lines are the analytical prediction of 
the scaling theory, without any adjustable parameters. (Only 
the lines for e/A — 15, 20, and 26 are shown for clarity.) The 
dashed curves are the golden- rule prediction (P), with a single 
adjustable parameter (the same for all curves, but the data 
for e/A = 15 was left out of the fit). 

The large- a; regime is described by the golden rule 
-^goidcn-ruie = S/T, according to which all basis states 
within the decay width F of a non-interacting state are 
equally mixed into the exact eigenstate. This complete 
mixing amounts to fully developed chaos. For our model 
the level spacing of the many-particle states in the j-th 
layer is (5 ^ j^^'^A/V{j) and the Breit-Wigner width is 
F ~ V^/A2 - f/^g-^A, which leads to Eq. (|). One 
notices in Fig. |l| that for the largest x the data points 
fall somewhat below the golden-rule prediction. This is 
due to the finite band width of the layer model. The 
IPR saturates at 3/'P{j) ^ when the decay width F be- 
comes comparable to the band width j^^^A. The corre- 
sponding upper bound on x for the validity of the golden 
rule is a; ^ j^^^lng. The finite band width of the layer 
model becomes less significant for large j, which is why 
the agreement with the golden rule improves with in- 
creasing j. 

The small- a; regime is described by the scaling func- 
tion F{x). The term of order a;" in the Taylor series (||) 
contains the n -I- 1-th order effective interaction V^^^ be- 
tween n + 2 particles and holes. A Fock state in the j-th 
layer contains about V7 excited particles and holes 
Because this is a large number for j ^ 1, the IPR fac- 
torizes into a product of independent contributions from 
2,3,4,... interacting particles, 

n=0 

A calculation of \V^fi \ leads to Eq. (|). The appearance 



of the modulus of the matrix element in Eq. (^) is easily 
understood for the case of only two unperturbed many- 
particle states interacting via the matrix element V^^ . 
The IPR changes by order unity if two Fock states come 
energetically within a separation \V°^\ of each other. 
The probability of such a near degeneracy is small like 
|]/°*^|/A. (There is no level repulsion for the many- 
particle solutions of the non- interacting Hamiltonian.) 
Because for weak interaction the IPR can change signif- 
icantly but only with a small probability, the IPR fluc- 
tuates strongly. Indeed, in our simulations much larger 
statistics was necessary in order to reach good accuracy 
in the small- a; regime. (The remaining statistical error in 
Fig. 1^ is smaller than the size of the markers.) 
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FIG. 2. Averages — In J, — In /, and In 1// as a function of 
g, rescaled in the same way as in Fig. 0, for e/A = 20. For 
small X, the three averages follow the scaling theory (Q) (solid 
lines). For large x the averages — InJ and Inl/J follow the 
golden rule (dashed line). 

In Fig. H we compare the logarithmi c ave rage In J with 
the two other averages In / and — In 1 //. Within the 
small- a; regime of validity of Eq. (||) the three averages 
are related by 

S7=2(2-\/2)lnT=-2(\/2-l)lnT7/. (7) 

These numerical coefficients do not depend on the num- 
ber of particles involved in the interaction (cf. the explicit 
calculation of / in Ref. [@). As one can see in Fig. ||, 
for e/A — 20, the relation (Q) agrees well with the sim- 
ulation. In the chaotic regime, for large x, Eq. is 
no longer valid. The average — Inl//, which is domi- 
nated by the majority of states having a large number of 
components, is close to In I at large x. The average In/ 
is dominated by rare states with an anomalously small 
number of components and falls below the two other av- 
erages. This indicates an asymmetric distribution of In/ 
in the chaotic regime for the layer model. 
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So far we have only addressed the question of the scal- 
ing variable that governs the transition to chaos. What 
remains is the question: How sharp is the transition? The 
singular threshold predicted in Ref. jlj] developes only 
in the thermodynamic limit and would be smoothed by 
finite-size effects in any simulation. The corresponding 
non-analyticity of In / is related to the high-order behav- 
ior of the series (^). Since our numerics allows to distin- 
guish only the first two coefficients cq and ci, it leaves 
open the question about the non-analyticity. Still, even 
if the series (|^) would be absolutely convergent, the re- 
sulting smooth function of the single variable x could not 
describe the IPR for large x because it is incompatible 
with the golden rule — j/^^ln/goidcn-ruic ^ x^^ In*?. This 
different scaling behavior for small and large values of x 
suggests that the peak observed in Fig. || would evolve 
into a singular threshold in the thermodynamic limit. 
The only way to maintain a smooth crossover would be 
to introduce a parametrically large interpolating region 
between the two different scaling regimes. We can not 
exclude this interpolating region on the basis of the nu- 
merical data, however, theoretically there is no indi- 
cation for such a region. 

In summary, by exact diagonalization of a model 
Hamiltonian we have presented evidence for an 
interaction- induced transition to chaos in a quantum dot. 
Upon inclusion of finite-size effects, a good agreement is 
obtained with the scaling theory of Ref. |14| , support- 
ing the assertion that x — (e/.gA)lng is the scaling pa- 
rameter for the transition. The different behavior of the 
scaling function for small and large x suggests that the 
transition would become a singular threshold in the ther- 
modynamic limit. 
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